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Spherical,  Cylindrical  and  One-Dimensional  Flows 
of  Compressible  Fluids 

Joseph  B.  Keller 


I . Introduction 

Relatively  fe\>f  boundary  value  problems  involving 
spherical  or  cylindrical  flows  of  compressible  fluids  have 
been  solved  exactly.  In  order  to  solve  more  problems  of  this 
type,  particularly  those  involving  variable  entropy,  we  have 
investigated  the  flow  differential  equations,  at  first 
without  regard  to  initial  or  boundary  conditions.  In  this 
way  a class  of  non- isentropic  solutions  of  the  differential 
equations,  depending  upon  an  arbitrary  function,  has  been 
obtained.  Then,  in  the  second  phase  of  the  investigation, 
the  arbitrary  function  is  adjusted  to  satisfy  particular 
initial  or  boundary  conditions.  In  this  way  the  free  expansion 
of  a sphere  of  gas  into  a vacuum  has  been  treated,  as  well  as 
the  propagation  of  finite  and  strong  shocks  in  variable  media. 
The  latter  treatment  includes  Frimakoff's  point-blast  solution 
as  a snecial  case. 

The  method  of  procedure  is  quite  simple.  In  section  II 
the  problem  is  formulated  in  Lagrangian  variables.  In 
section  III  a class  of  solutions  is  obtained  by  the  method  of 
separation  of  variables.  In  section  IV  these  solutions  are 
specialized  to  isentropic  flow.  In  section  V the  solutions 
are  applied  to  the  do  termination  of  strong  shock  waves 
propagating  in  variable  media.  In  section  VI  finite  shocks 
in  variable  media  are  considered. 

II.  Formulation 

VJe  consider  the  mo  cion  of  an  inviscid,  non-heat- 
conduct-ing  fluid,  obeying  the  polytropic  equation  of  state. 

The  one,  two  and  three  dimensional  cases  will  be  treated 
together.  In  the  three  dimensional  case  y(h,t)  represents 
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the  radius  at  time  t of  the  particle  with  the  Lagrangian 
coordinate  h,  which  is  defined  by  the  equation 


(i) 


-I 


rn_1p(r,t)  dr 


n = 1,2, 3- 


y( o,t) 


In  the  above  equation  p(r,t)  is  the  density  at  time  t and 
radius  r,  and  n is  the  dimension,  which  is  3 in  the  spherical 
case.  In  the  two  dimensional  case  y represents  radial 
distance  from  an  axis,  and  in  one  dimension  y is  a cartesian 
coordinate . 

It  is  further  assumed  that  all  flow  variables  depend 
upon  y and  t only  and  that  flow  occurs  only  in  the  y direction. 
This  is  the  assumption  of  spherical,  cylindrical  or  planar 
symmetry . 

From  the  definition  of  y,  the  particle  velocity  u is 
given  by 


(2) 


u = y. 


Similarly  from  (1) 
given  by 

(3) 


the  density  p or  specific  volume 


n-1 

y y 


h 


If  are 


Because  of  our  assumptions  concerning  inviscidity  and 
non-conduction,  the  entropy  s of  a particle  is  independent 
of  time  (at  least  between  successive  shocks).  Thus  we  have 


( 4-)  s = s(h) 

The  function  s(h)  is  given  by  initial  data  or  by  shock 
conditions,  and  is  assumed  to  be  known. 

The  pressure  p is  given  by  the  equation  of  state 

(5)  P = P(p, s)  = g(V,s) 


For  a polytropic  gas  or-  liquid 
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(6)  g(t?,s)  = g0  + A(s)2fY 

The  function  A(s),  the  adiabatic  exponent  y>  and  the  internal 

pressure  g are  assumed  to  be  known, 
o 

In  terms  of  the  above  defined  quantities,  the  equation 
of  motion  is 

(7)  ytfc  = -yn“1fgu(yn‘1yh)h  + sssh3  . 

For  a polytropic  gas  or  liquid  this  becomes,  using  (6), 


(8) 


y 


tt 


YA(  s 


^ ( 


n-1  v - y- 1 , n-1  v 

t yh)  1 (y  yh) 


n-1 


, n-1  \ — Y 

(y  yO 


*h 


n-1 


h 


Equation  (8)  is  a second  order  partial  differential  equation 
for  y(h,t).  The  coefficient  A is  assumed  to  be  a known 
function  of  s,  and  by  ([;.),  of  h. 

The  problem  wc  consider  is  that  of  finding  solutions  of 
(8).  From  a particular  solution,  the  flow  variables  can  be 
found  by  using  (2),  (3),  (l\.)  and  (5)* 


III.  Product  Solutions 

Let  us  seek  product  solutions  of  (8)  of  the  form 
(9)  y(h,t)  = f(h)j(t) 

Inserting  (9)  into  (3),  and  separating  variables  we  obtain 


(10) 

(11) 


-A!(fn-1f'>-q,fn-2  -A,;fn-1f'rY  fn"2  = x 


In  these  equations  A is  an  arbitrary  separation  parameter  and 
primes  represent  derivatives  with  respect  to  t in  (10)  and 
with  respect  to  h in  (11). 

To  solve  (10)  we  multiply  by  j ’ and  integrate,  obtaining 


(12) 


1/  x 2 A ,n(l-v)  _ 1 

J } ' nTT-77  j " 2 


(Y  t 1) 


k 


(1.3)  £■(  1 ) - X lo  : j - a (y  - 1)  * 

Hero  a is  an  integration  constant.  Thus,  unless  j is  constant, 
which  is  possible  only  if  X = 0 or  j ~ 0,  we  have  the  solution 

04)  J [hTT^TT  jn(1'Y>  +a  1/2dJ  = t (y/1) 

(15)  f (2X  log  ) + aj'1//2  dj  = t (y  = 1) 

The  integration  constant  in  (1)+),  (--5)  has  been  sot  equal  to 
zoi-c  wither • ’’oyp  of  generality  since  it  merely  shifts  the 
origin  of  time. 

To  solve  (11)  we  differentiate  out  and  obtain 

(16)  YAlfiifn"1  + ( n-1 ) fn-1(  f 1 )2)(fn"1f ' J'Y-1 

- fn“2(f11_1f ' )~y  A*  = X 

We  now  consider  the  inverse  function  h = h(f)  and  denote 
h'(f)  by  q(f).  Then  (16)  becomes 

(17)  YA[-qV3fn_1  + (n-l)fn"2q"2][fn“1q"1]"Y"1  fn~2 

+ fn_2r J_Y  A ' ( h ) = X 

We  shall  first  treat  the  case  y ^ 1 by  introducing  z(f)  and 
3(f)  defined  by 

(18)  z = qY-i  , B(f)  = A[h( f ) j 

With  these  definitions,  (17)  becomes  (prime  denoting 
differentiation  with  respect  to  f) 

(10)  z'  + z[-(n-l)(Y-l)f'1  + I~  (logB)'] 

M Y-l)  f(n-l)(  y-D+1  _ 

yb 


+ 


0 


and  the 


An  integrating  factor  of  (19)  is  f Y 

solution  is  (with  G a constant) 


llX 


f 


(20)  a = f(n-i)(r-i)B  y tG  . C fB"i/T  an 

Thus  we  have  q from  (18)  and  (20),  and  finally  since  h’  = q, 

r 


(2i)  h = J fH-ig-i/r  g - M.rrAl  f 


fB 


-1/y 


df 


Y-l 


df  + E 


Here  5 is  an  integration  constant. 

Ec”--<tion  (21)  gives  f(h)  implicitly,  and  thus  provides 
a solution  of  (11)  for  y / The  solution  may  be  written 

more  simrly  by  defining  P(f)  by 


(22) 


’(f)  = G - ALllii  J 


fe 


-i/y 


df 


Y-l 


If  \ f 0 this  can  be  solved  for  B(f)  and  yields 


(23) 


tdJ  f* 


(-xf)Y(  f’)~Y  p 


The  solution  for  the  flow  variables  can  now  be  computed  from 
(2) -(5),  (9),  (21),  (22).  We  have  for  y / 1 and  \ / 0, 
remembering  that  f = yj~”  from  (9), 


(21*) 

(25) 

(26) 


u(y,t)  = yj'j"1 

^(y,t)  = -Xyjn“1/F’ (yj*1) 

p(y»t)  = gQ  + j"nY  ^‘(yj"1) 


Equations  (2k) -(26)  give  the  flow  quantities.  In  fact,  these 
expressions  yield  a solution  of  the  Eulerian  equations  of 
motion  for  an  arbitrary  function  p provided  that  j(t)  is  given 
by  (ll|-).  This  is  the  first  main  result  of  this  paper.  It 
is  to  be  noted  that  in  these  solutions,  u is  proportional  to  y. 

In  the  excluded  case  y / k ^ = 0 we  have  instead  of 
( 2k) -(26)  from  (2)-(5),  (9),  (21) 
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(27) 

(28) 

(29) 


u(y,t)  = yj’j 


-1 


= yt 


-l 


r(y,t)  = jnB1/Y(yj"1)G1"Y  = tnb(yt-1) 
p(y,t)  = g„  ♦ rnr  a*  = . 


The  new  arbitrary  function  b and  arbitrary  constant  S.  have 
been  introduced  and  the  expressions  simplified  by  using  ( II4.) , 
which  gives 


(30) 


j(t)  = a1//Pt 


The  corresponding  solutions  when  y ~ 1 are  obtained  from 
(17)  by  introducing  3(f)  as  in  (18),  after  which  (1?)  becomes 


I -1  , , . _-l  , 3 . Xf  _ , 
q q -(n-l)f  + ~B~  + ~b  ° 


(3D 

The  solution  of  (31)  is 

rf 

(32)  q(f)  = fn_13"1(f)  exp  J -Xf3~1(f)df  . 

r 


Thu  s 

(33)  h(f) 


fn_13’1(f)  exp  J -Xf3_i(f)df 


df  + E 


Equation  (33)  yields  the  solution  implicitly.  Again 
we  define 

( 3U-)  2(f)  = exp 

Then  if  X f 0 we  have 
( 33)  3(f)  = -\f F(F ' ) 

The  solution  for  the  flow  variables,  if  y = 1 and  X ^ 0,  is 
again  given  by  ( 2Lp)  - ( 26 ) with  y = 1 and  j(t)  given  by  (19). 
Similarly  for  y = 1 and  X = 0 the  solution  is  given  by 
(27) -(29)  with  y = 1,  G = 1 and  j(t)  given  by  (30). 


7. 


Equations  (2[l) -(26)  and  (27) -(29)  represent  non- 
Isentropic  solutions  of  the  flow  equations  depending  upon  an 
arbitrary  function.  In  order  to  construct  these  solutions 
explicitly  one  need  merely  evaluate  the  integrals  in  (Ik) 
or  (l£).  Some  special  cases  of  these  integrals  yield  results 
which  are  listed  below  for  later  use. 


36)  = 


if  y 1,  a = 0,  yi  _i 


(37) 

±x1/2t 

l(  t)  - e 

if  y/  1,  a=  0, 

n ( y- 1 ) _ -i 

2 

(33) 

j(  t)  = ^ t2+  at  + fi 

if  n(l- 

•y)=  i. 

a,  p arbitrary. 

IV. 

Isentropic  Case 

The  above  solution 

can  be  specialized  to  the  isentrcpic 

case 

by  setting  A = 3 = 

constant.  Then  from  (22) 

and  ( 3k)  we 

obtain 

(39) 

F(f ) = [ G - 

1 2TB1/Y 

p]  Y-l 

n y ¥ i 

(ko) 

_2 

F(  f ) = exp 

Y = 1 

• 

Thus 

the  solutions  given 

by  (2k) -(26)  become,  for 

Y / 1 and 

X = i 

(4D 

u(y,t)  = yj' j"1 

r 

■1 

(k2) 

V(y,t)  = j^b1^  [ g 

- ihjU  yh-2]B 

2yB  /Y  J 

■1 

(43) 

p(y»t)  = gc  + rnY  | 

'o  - ^ y2J'21 

2YBX/Y 

V 

Y-l 

• 

In  these  solutions  3 and  G are  constants  -and  j is  given  by  (Ik). 
For  X = 0,  and  all  y The  solution  (27)-(29)  applies  with  b 
constant , 
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For  y = 1 and  all  (41)  is  unchanged,  but  the  other 
equations  become 

-2  2 

(44)  tf(y,t)  = jnB  exp 

-2  2 

(45)  P(y,t)  = g0+  j"n  exp 

Here  3 is  a constant  and  j is  given  by  (15)* 

As  an  example  of  these  isentropic  solutions,  let  j(t) 
be  given  by  (36)  and  assume  G = 0 in  (42),  (43)*  Then  (4*1)- 
(43)  yield  the  power  solutions  (if  y / 1,  — - '-Cf----  / -1) 

(46)  u(v,t) 

(47)  «(y,t) 

(48)  p(y,t) 


V • Application:  Strong  Sho cks  In  Variable  Media 

Let  us  suppose  that  a shock  given  by  the  equation 
y = R(t)  moves  into  a variable  mediun  of  density  pQ(y), 


cond  it ions  s 


(49) 


SL.  = 

P 


(.50)  u = 


(5D 


R = 


,(y)  and 

velocity  zero.  The 

pres 

' p just- 

behind  the  shock  i 

are 

rel, 

.ng  quantities  in  front  of 

it 

by 

These 

conditions  are 

( Y+1)P  + 
( Y-1)P + 

(Y-1)po  Y+l 

TFTrr  * trr 

(Y  i 

S(P  -p0) 

frv 

{,2p0[(  Y+l)p+  ( y-Dpo]}1//2 

L(7 

f ( r+Dp 

+ (Y-l)po1l/2^  r 

( Y+Dp] 

2Po  J L 

2 

Po  J 

11/2 


1/2 
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The  second  expressions  on  the  right  apply  to  a strong  shock, 
for  which  p » pQ. 

We  shall  now  assume  that  the  flow  behind  the  shock  is 
a "product"  solution  given  by  (24) -(26),  and  that  the  shock 
is  strong.  We  wish  to  determine  the  functions  F(yj-1),  R(t), 
po(y)  and  pQ(y)  for  which  such  a solution  is  possible.  To 
this  end,  we  insert  ( 2I4- ) — ( 26 ) into  (49)  -(51)  and  obtain 


(52) 


FllSilil  = V+l 

-XRj^"1p0(R)  ^ 


(53) 


Pp*  - 


2j“nY?(Rj' 


Y 


(Y+l)fg0+  j’nYF(Rj_1)} 


(54) 


p = 


1/2 


Prom  (53)  and  (54)  wo  have 


(55) 

Thus 

(56) 


R _ (v+l)  j 

R ’ 2 ~ 


R(t) 


xii 

R_,[  j(  t)  ] 


Inserting  (56)  into  (53)  and 

making 

/ use  of  (12) 

> 

sat  i 

sfied  by  j(t),  we  find 

/ •Y— — —6 — A 

(57) 

F!x!  - Po  (x^  ,l-r) 

■ i-x 

-2ny 

( 2n+2 ) y-2 

2ny 

IRJ’1 

y- 1 
X 1 

Thus  D is  determined,  F is  related  to  p , and  (53),  (5U-)  are 
satisfied,  although  Rq,  P,  and  >.  are  still  arbitrary.  We  now 
insert  (56),  (57)  into  (52)  to  determine  F.  We  obtain,  if 


o 


10 


(58)  F(x)  = 


F x 
o 


n 


2XR2n  ? 1-1/2 
o , -,„2n 

Tn^rT  M 


Froni  (57)  and  (58)  we  then  find 


2F  R? 


£21 

y-1  r 


o o 


(59)  P0(R)  = y+! 


2XR 


2n  Jm  2n(  y-1)  ~) 

t + dr^-  r Y 1 

n(l-y)  o 


-3/2 


In  (58)  and  (59)  FQ  is  an  arbitrary  constant.  The  corres- 
ponding results  with  gQ  / 0 are  somewhat  more  complicated* 

We  have  thus  obtained  a solution  with  a strong  shock 
moving  into  a variable  medium  at  rest,  with  density  given  by 
(59)-  The  constants  Fq,  Rp,  \ and  D are  arbitrary  in  this 
equation,  but  only  two  essential  combinations  of  these 
constants  occur.  The  shock  curve  is  given  by  (58),  and  j(t) 
by  ( lip)  . The  flow  is  given  by  ( 2ip)  — ( 26 ) with  F(x)  given  by 
(58).  The  flow  might  be  produced  by  a piston  following  one 
of  the  particle  paths 


(80)  Y(t)  = Y0.j(t) 

The  solution  was  deduced  for  y ^ tl  and  X / 0. 

As  an  example  of  these  solutions,  lot  us  suppose  that 
P0(R)  = constant.  From  (59)  we  find  that  this  requires 

(61)  »“<>  . T-^|  • 

The  last  condition  can  be  fulfilled  only  for  n = 3,  in  which 
case  y = 7.  (n  = 1 leads  to  y = -1,  which  was  excluded  in 
the  derivation,  and  n = 2 yields  no  value  of  y. ) In  this 
case  we  have  from  (!/.(.),  (58) 
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Tho  solution  computed  by  using  (62)  in  (24) -(26)  (with  X < 0) 
is  exactly  the  point  blast  solution  of  Taylor’s  type  first 
found  by  H.  Primakoff, 


VI . Application:  Finite  Shocks  In  Variable  Media 

We  may  apply  the  above  method  to  determine  finite 
shocks  in  variable  media.  Then  we  must  satisfy  the  exact 
shock  conditions  (49) -(5l)  rather  than  the  strong  shock 
conditions.  Inserting  ( 2ip)  — ( 26 ) into  (49)-(5l)  yields  with 

g =0 

- o 

(rH)rnYF+(Y-i)Po 

— = , x = p j 

p0  (Y+1)P0+  (y-1)  3 YP 

i'nrF  -p 

(64)  Rj  j 1 = — 2 

j^[  ( y+1  ) j”nYP  + ( y-1 ) P0 1 1/2 


(63) 


-XRj 


L.l*l 

, .n-1 


(65)  R = — - — [(y+1)  j'nrF^  (Y-l)po}1//2 

Equations  (63) -(65)  are  a set  of  two  first  order 
ordinary  differential  equations  and  one  algebraic  equation 
for  the  determination  of  tho  four  functions  F(x) , R(t),  PQ(R) 
and  pq(R).  Equation  ( 14)  gives  j(t).  This  system  is  evidently 
underdetermined  and  will  have  infinitely  many  solutions.  In 
the  strong  shock  case,  however,  p did  not  occur  and  thus  the 
system  was  determined. 

To  solve  the  above  system  wo  could  impose  some  relation 
between  pQ  and  pQ  and  eliminate  both  these  functions  by 
means  of  that  relation  and  (64).  A pair  of  simultaneous 
first  order  equations  for  F(x)  and  R(t)  would  result.  However 
these  equations  can  be  treated  separately  since,  from  (64) 
and  (65)  we  have  by  eliminating  F, 


12. 


(66)  r = ± [iX+li!(Rj'3-l)2+YpoP-l]1/2  . 


-1 


This  is  an  equation  for  R(t)  if  the  ratio  P0P0  is  given  as 
a function  of  R.  (Equation  (66)  also  holds  when  gQ  ^ 0.) 
After  solving  this,  (63)  can  be  solved  for  P,  and  then  pQ 
and  p can  be  obtained. 

We  will  now  investigate  the  shock  curve  R(t)  given  by 
(66),  in  the  special  case  in  which  j(t)  is  given  by  (36)  with 
C = 0 and  YP0P~^  = cq  ^-s  constant.  The  latter  assumption 
means  the  temperature  ahead  of  the  shock  is  constant.  Making 
use  of  (36),  (66)  becomes 


(67) 


A _ Y+ 1 - - 1 

R - ?[n( y-IT+21  rtt 


(v+l)‘ 


Ii4-[n(  y-1)+2] 


■1 


2 (Hf1)2  + C2 


1/2 


If  we  now  introduce  TT(t)  = Rt“  in  (67)  we  obtain  an 
equation  for  U in  which  the  variables  separate.  One  solution 
is  the  constant  solution 


(68) 


TT  = + 


+ „ ,,  Y+l  >-1/2 
1 


In  this  case,  which  is  physically  possible  only  when 

—7 — t , -v  < 1,  the  shock  curve  is  the  straight  line 
n(  v- 1 ) 


(69) 


R = nt  = ± ‘ 


When  W does  not  have  the  value  given  by  (68),  the 
equation  for  TJ  yields 


(70) 


U(t) 

JIt  {‘srPT 


- 1)11 


±f hlih  ,,2  + ^l172]-1 


u [2n(Y~l>+43‘ 


°J  ) 


dU  = log 


Integrating  (70)  yields 
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(71)  bt  = i a + (a-l)cos  9 


(sin  fj^'1  (cos  |^a 


_1 

-a+1 


-1  aU 


Here  9 = cot  ^ , a = -T-7-'i— 1 and  b is  a constant.  Thus  we 
c 2n(  y-D+4 

have 


(72) 


R = T.Tt  = 


bt  cot  9 


^ince  bt  is  given  in  terms  of  0 by  (71)>  R(t)  is  given 

fit) 

parametrically  in  terms  of  9.  A graph  of  — R versus  bt  is 

co 

given  in  Figure  1 for  n = 3 and  y = 1»4»  which  represents  a 
spherical  shock  in  air.  When  the  minus  sign  is  chosen  in 
(67) -(71)  the  same  curves  are  obtained  as  with  the  plus  sign, 
provided  t is  replaced  by  minus  t. 

For  small  values  of  t we  have 

(73)  R » Kt2a 


where  K is  a constant.  For  y = 7 and  n = 3 we  find  a = -r-  and 
2 /[3  b 

thus  R ~ t ' , which  is  the  behavior  of  Primakoff  ’ s point 

blast  solution  for  small  t.  For  large  t,  R behaves  linearly 

in  t,  as  in  {69). 

To  determine  pn  and  F,  we  have  from  (65)  and  the  condition 


YP0P0  = co’ 


(7U) 


2Pr!T  - ( y-1)  p 


h+l)j"nY  ( y+l)  c' 


° R2inY  . Hi  ,ny 


y+l 


Thus,  if  dot  denotes  time  derivative  and  prime  denotes 
derivative  of  a function  with  respect  to  its  argument,  we  have 

(75)  F ,-4r_^(p;s3,nY+  KjnY+ 

( Y+Dc 


■ ^l(Po^'nY+  nYFc 


.ny-1 . 1 


T , * - 1 _ p t 

= F { R j 1 -Rj  j ) 


4 sing  ( Tn ) , (75)  in  (63)  yields  the  following  equation  for  p 
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(76) 


L(  y+1  ) c 


<P^3JnY*  + ) 


o 


- Sr(P^.1nY^YP0.in^1j')l  ["-XR3n"1YP0c;2(Hj'1-  Rj2j') 


1-1 


J 


2Yc-o2R2 


vfl+  (y-i.).2yR2  (y-D.2 

( Y+l) co  Y+1  - 


-1 


Upon  introducing  = p^R,  the  above  equation  can  be  integrated 
yielding  p^[R( t) ] as  a function  of  t.  To  obtain  p (R),  we 
must  insert  t = t(R)  from  (71) » (?2)  into  this  result.  Then 
when  p [R(t)j  is  known,  F can  be  found  from  (7 lx).  These 
calculations  will  not  be  carried  out  here. 


V T I . Free  Expansion  and  Other  Problems 

In  this  section  some  flows  described  by  the  preceding 

"product"  solutions  will  be  discussed.  First,  let  us  consider 

the  flow  given  by  ( Lpl ) — ( Ip 3 ) with  gQ  = 0.  This  describes  the 

isentropic  flow  of  a pras.  If  j'(0)  = 0,  it  starts  from  rest. 

If  the  constant  G is  chosen  so  that  p = 0 for  some  value  of 

yj-^,  then  p remains  zero  for  this  value  of  yj' ^ , which 

therefore  represents  a "free"  surface  of  the  gas.  Thus  the 

flow  represents  the  free  expansion  of  an  isentropic  gas  into 

vacuum.  The  equation  of  the  free  surface  is  y = j(t), 

where  v represents  the  location  of  the  free  surface  at  t = 0. 
o 

In  a similar  way  (2)4.) -(26)  may  represent  the  free  expansion 
into  vacuum  of  a gas  with  initially  variable  entropy,  provided 
gQ  = 0 and  F(yj-'  ) vanishes  for  one  value  of  its  argument. 

For  certain  choices  of  j(t),  this  solution  can  also  represent 
the  collapse  or  expansion  of  a hole,  or  the  initial 
contraction  of  a sphere  of  gas  toward  the  origin  due  to  an 
initial  inward  motion,  followed  by  momentary  rest  and 
subsequent  expansion.  When  gQ  is  not  zero,  similar  motions 
of  a liquid  are  described. 
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VIII.  Other  Solutions 

In  the  preceding  sections  product  solutions  of  (3) 
were  considered.  In  this  section  two  other  kinds  of  solutions 
will  be  examined.  First,  let  us  assume  that 

(7?)  y = f ( at  + ph)  . 

Here  c and  6 are  constants.  Me  find  for  the  flow  quantities 
(73)  u = = cf  = cf [ f_1(y) ] 


(79)  ^ = p-1  - 6e-1yn_1u 

(30)  P - ApY  = A(ap'1y^“nu~1)r 


Me  see  that  the  flow  quantities  are  independent  of  t,  or  In 
other  words,  that  a solution  of  the  form  (77)  yields  steady 
flows. 

Inserting  (77)  Into  (3)  yields 


c f = yA( 


nr11 


in. 1 • * 

(?fa  ""’f ) 


pr1'1  + ,\h(  Qfn_1f )' 


Yfn 


-1 


Since  f is  a function  of  ct+ph  and  t doesn't  appear  in  (3l), 
h cannot  sop-ear  either.  Thus  A is  constant  and  = 0.  Now 
if  x denotes  the  argument  of  f,  and  if  we  consider  the  inverse 
function  x(f),  we  obtain  from  (ril)  the  following  equation 
for  q(  f ) = x(f) 


(32)  -c2qq-3  = YA(pfa  Xq  1}”Y“1(-6fn  ^q”3  + B(n-l)fn“£:q“*)prn“-L. 


-n-2  -2A„,.n-l 


This  equation  is  exact  and  yields  upon  integration 


(8l)  bf 


( 1~y) ( n-1)  Y-1  1 -2  _ 


c 


2 « 


yApJ 


.1—. 


c-(y-l) 


(3M 


-2 


31;)  -j-  q + A log  q = (n-1)  If  +D 


T ¥ 1 
Y = 1 


H-  r-e  v is  a constant- 


16. 


PrCi"  ( 7-2 ) and  the  Tacts  that  f = 


-1 


f = v.  we  have 


(S5) 


u 


qlyT  * 


Thus  the  steady  solutions  are  obtained  from  the  algebraic 
equation  (8l)  or  {3!;.}. 

A more  interesting  type  of  solution  is  given  by 


(86)  y = taf(K(h)t) 

Here  a is  a constant,  and  f,K  are  functions  to  bo  determined. 
These  solutions  aKr  like  the  type  considered  by  Guderly, 
Taylor,  von  Pe unarm  and  Calkin  and  Courant  and  Priedrlchs. 
Inserting  (36)  Into  (5)  yields 


{37)  (a-l)ata"2f + 2aKta"if + K 2V 


- yA  1 1 


ftan+1 ^-1  fp,  Y~V  ( n-1 ) f n-1  fcan+l 


frn-i 

L 


nc  + 


rw  1 * “P  mm 

f-  Vtf  + (n-1 


+ 


In  order  that  h occur  in  this  equation  only  in  the  combination 
x — I-  ( h ) t , we  find  the  following:  expressions  for  K ( h ) and  A(h) 


{53}  K = ^3h+C , A = DKy-1Kg+1 

Here  c - an{ 1-y  s -y-2n+2,  D is  a constant  and  A is  constant  if 
e(y-l}+e+l  - C . 


Tinder  these  conditions,  (8?)  becomes 


a-  ■ 


(39)  (a-l)af + 2axf + x“f  = 


xcf{n-l)(l-Y>^-rD 


e + 


xf  (n-l)xf 

T 


£ { y— 1)  + r + 


1?. 


To  reduce  (39)  to  a first  order  equation,  v:e  follow 


-1 


von  Neumann  and  Calkin,  by  letting  6 = ( a+y)  [2  f n(  y-1)  ] and 


(90) 


x = es  , f = e^sF 


Then  (89)  becomes 


(91)  (a-l)aF+  2a(^  + 6)F  + (£  + S-l ) (£  + 5)F 


as 


ds 


= P 


( n - 1 ) ( 1 - y ) 


d 1 "Y 

{—  + 6)F  D 
. ds  J 


ye 


(-£_+  6 ) f 
d s 


Y(n-l)  (-?-+•  6)F 

+ 3^ + ey-e  + c+1 


If  we  now  let  q = {^  + 6)F  we  obtain  a first  order 
equation  for  q in  terms  of  F. 


(92)  |l-  y^F 


( n-1 ) ( l-y_)  -y-1 


( q - 6F)  q 


= F(n_1' ’ 1_Y)q"YD[  Yt  + Y(n-l)qF"1+ey+a+l-e]-(a-l)aF-2aq 


Special  cases  of  this  equation  have  been  studied  by 
von  Neumann  and  Calkin. 


Figure  1 
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Because  of  our  ] 
YOUR  PURPOSE 
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